Magnetic properties of BaCdVO(P04)2: a strongly frustrated 
spin-1/2 square lattice close to the quantum critical regime 
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We report magnetization and specific heat measurements on polycrystalline samples of 
BaCdVO(P04)2 and show that this compound is a S = 1/2 frustrated square lattice with ferro- 
magnetic nearest-neighbor (Ji) and antiferromagnetic next-nearest-neighbor (J2) interactions. The 
coupling constants Ji ~ —3.6 K and J2 — 3.2 K are determined from a fitting of the susceptibil- 
ity data and confirmed by an analysis of the saturation field (fioH s — 4.2 T), the specific heat, 
and the magnetic entropy. BaCdVO(PC>4)2 undergoes magnetic ordering at about 1 K, likely to- 
wards a columnar antiferromagnetic state. We find that BaCdVO(P04)2 with the frustration ratio 
q = J2/J1 — —0.9 is closer to a critical (quantum spin liquid) region of the frustrated square lattice 
than any of the previously reported compounds. Positive curvature of the magnetization curve is 
observed in agreement with recent theoretical predictions for high-field properties of the frustrated 
square lattice close to the critical regime. 

PACS numbers: 75.50.-y, 75.40.Cx, 75.30.Et, 75.10.Jm 
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I. INTRODUCTION 

Low-dimensional spin systems are one of the actively 
studied subjects in solid state physics due to the pos- 
sibility to observe numerous quantum phenomena and 
to interpret these phenomena within relatively simple 
models (e.g., Ising or Heisenberg models for different 
lattice types). An interesting phenomenon in the spin 
physics is the formation of a spin liquid - a strongly 
correlated ground state lacking long-range magnetic or- 
der. This ground state is usually related to the electronic 
mechanism of superconductivity suggested for high-T c 
cupratesji^ Spin liquids originate from quantum fluctua- 
tions that are particularly strong in systems with reduced 
dimensionality and low spin value. The fluctuations can 
be further enhanced by introducing magnetic frustration 
which impedes long-range ordering of the system. 

The spin-1/2 frustrated square lattice (FSL) is one of 
the simplest models giving rise to a spin liquid ground 
state. In this model (also known as the J\ — J2 model), 
magnetic moments on a square-lattice are subjected to 
nearest-neighbor interaction J\ along the side of the 
square and next-nearest-neighbor interaction J2 along 
the diagonal of the square. The FSL model is de- 
scribed by the frustration ratio a = J2/J1 or, alterna- 
tively by the frustration angle ip =tan (J2/J1). Defin- 
ing the thermod ynamic energy scale of exchange cou- 
plings as J c — \J J \ + J 2 i on e obtains J\ = J c cos ip and 
<h = J c sin ip (see Fig. Q] and Ref. y|). 

Extensive theoretical research on the FSL model has 
been done in the past. Initially, the studies were fo- 
cused on the AFM region (Ji, J 2 > 0)££&£ while the 
general case (arbitrary signs for J\ and J2) was only 
considered quite recently 10 1 1 1 The phase diagram 
of the model (Fig. [1]) includes three regions with dif- 
ferent ordered phases: ferromagnet [FM, wave vector 
Q = (0,0)], Neel antiferromagnet [NAF, Q = (tt : tt)}, 



and columnar antiferromagnetic [CAF, Q = (tt, 0) or 
(0, 7t)]. Classically, first-order phase transitions should 
occur at the NAF-CAF {a = 0.5) and the CAF-FM 
(q = —0.5) boundaries. However, quantum fluctua- 
tions destroy long-range ordering, leading to the forma- 
tion of critical regions with disordered ground states. In 
general, the ground state in the critical regions is re- 
ferred to a quantum spin liquid (QSL) regime, but the 
particular nature of the spin liquid phases is still un- 
der discussion. A gapless nematic state is suggested 
for a ~ —0.5, 8 while different dimer phases (including 
resonating-valence-bond-type ones) are claimed to exist 
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FIG. 1: (Color online) Phase diagram of the FSL model.— 
Solid filling indicates the regions of long-range magnetic or- 
dering, while hatched filling denotes the critical regions. Po- 
sitions of BaCdVO(P04)2 and some of the previously inves- 
tigated compounds are shown (see text for references). 



for a close to 0.5. 5 Note also that the boundaries of the 
critical regions are not known exactly: for example, an 
a range from 0.34 to 0.60 is reported in Ref. H for the 
spin- liquid regime, while a wider region (0.24 — 0.83) is 
suggested in Ref. 

Despite numerous theoretical investigations, experi- 
mental realizations of the J\ — Ji model are scarce. 
Layered vanadium oxides Li2VOX04 (X = Si, Ge) and 
VOM0O4 were the first examples of the FSL systems. 
Initially, these materials were ascribed to a ~ 1 region 
of the phase diagram, and frustration-driven structural 
distortions were conjectured on the basis of nuclear mag- 
netic resonance (NMR) data . 13 ' 14 ' 15 However, later stud- 
ies revealed the lack of the structural distortions and es- 
tablished a different scenario for all three systems. Thus, 
Li2VOX04 compounds fall into the CAF region with 
a > l ) 16 i 17 i 18 while VOM0O4 lies in the NAF region 
with a ~ 0.2«±2, 

Complex vanadium phosphates AA'VO(P04)2 present 
another realization of the FSL model with ferromagnetic 
Ji and antiferromagnetic J2 as will be discussed in detail 
below. We should also mention two very recent proposi- 
tions for the FSL systems. (CuX)LaNb 2 7 (X = CI, Br) 
perovskite-type compounds were claimed to realize the 
J\ — J2 model with ferromagnetic J\ and antiferromag- 
netic J21 20 ' 21 However, the estimates of J's from different 
experimental techniques are contradictory, and the valid- 
ity of the FSL model for these systems is still unclear^ 2 - 
The layered perovskite PbV03 also reveals an interesting 
square lattice system with a ~ 0.3, i.e., quite close to the 
antiferromagnetic (AFM) critical region . 23 ! 24 However, 
the complicated preparation procedure strongly hampers 
detailed investigation of this compound. Thus, little ex- 
perimental information about the critical regions of the 
FSL model is available, and the search for new FSL sys- 
tems is still challenging. 

Complex vanadium phosphates AA'VO(PC>4)2 have 
layered crystal structures (Fig. [2]) with square lattice-like 
arrangement of V+ 4 (S = 1/2) cations. [VOPO4] lay- 
ers are formed by VO5 square pyramids linked via PO4 
tetrahedra. The tetrahedra allow for superexchange in- 
teractions both along the side and along the diagonal of 
the square, hence the FSL-like spin system is formed. 
Metal cations and additional, isolated PO4 tetrahedra 
are located between the layers. The layers are flexi- 
ble towards buckling, therefore metal-oxygen distances 
can be tuned, and different metal cations can be accom- 
modated within the structure. Magnetic properties of 
the compounds with AA' = Pb2, SrZn, and BaZn have 
been recently investigated by means of thermodynamic 
measurement a 25 ! 26 and neutron scattering. 27 ' 28 We found 
that all the compounds fall to the CAF region of the FSL 
phase diagram. They reveal ferromagnetic J\ and anti- 
ferromagnetic J2 with a varying from —1.8 (AA' = Pb2, 
BaZn) to -1.1 (AA' = SrZn). 

An advantage of studying vanadium phosphates is the 
low energy scale for the exchange couplings (below 10 K) 
that allows for high-field experiments providing addi- 
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FIG. 2: (Color online) Crystal structure of BaCdVO(P0 4 ) 2 : 
single [VOPO4] layer (upper panel) and stacking of the layers 
(bottom panel). Arrows indicate superexchange interactions 
Ji (along the side of the square) and J2 (along the diagonal 
of the square) . Larger and smaller spheres denote Ba and Cd 
cations, respectively. 



tional information about the properties of the system^ 
Yet there is also one complication. The crystal struc- 
tures of AA'VO(PC>4)2 do not have tetragonal symme- 
try; therefore, vanadium atoms do not form a regular 
square lattice. At first glance, the distortion of the square 
lattice is negligible. However, even a very slight alter- 
ation of the structure can lead to drastic changes in the 
exchange couplings as shown recently for Ag2VOP2C>7 
(Ref. [23). In case of the AA'VO(P04)2 compounds, we 
studied this issue in detail using band structure calcula- 
tions and found considerable deviations from the square 
lattice model for some of the systems^ Nevertheless, the 
deviation is really negligible for one of the compounds - 
BaCdVO(P04)2- Below, we present magnetic properties 
of this compound and interpret the results within the 
FSL model. A detailed examination of the appropriate 
spin models for other AA'VO(P04)2 compounds will be 
published elsewhere j 3 ^ 

The crystal structure of BaCdVO(PC>4)2 has been re- 
ported by Meyer et aL^L but magnetic properties of this 
compound were not investigated. In our work, we use a 
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set of different thermodynamic measurements to study 
magnetic interactions in BaCdVO(P04)2- We show that 
BaCdVO(P04)2 is a FSL system with the frustration ra- 
tio a rs —0.9, i.e., lying very close the critical region of 
the FSL. 

The paper is organized as follows: Sec.|TT]dcals with the 
experimental details. In Sec. IIII| we present our results 
on BaCd(VO)(PC>4)2. Section ITVl contains a detailed dis- 
cussion of the results in the light of the J\ — J2 model 
followed by our conclusions. 



II. EXPERIMENTAL DETAILS 

Polycrystalline samples of BaCd(VO)(P04)2 were pre- 
pared by solid state reaction technique using BaC03, 
CdO, V 2 3 , V 2 5 , and (NH 4 ) 2 HP0 4 as starting mate- 
rials (all the chemicals had at least 99.9% purity grade). 
The process involved two steps. First, the intermedi- 
ate compound BaCdP2C>7 was prepared by firing the 
stoichiometric mixture of BaCC>3 and (NTL^HPCU at 
950 °C in air for 48 h with one intermediate grinding. In 
the second step, stoichiometric amounts of BaCdP2C>7, 
V2O3, and V2O5 were grinded, pelletized and annealed 
in dynamic vacuum (10 -5 mbar) or evacuated and sealed 
quartz tube (10" 2 mbar) at 800 °C for 30 h. 

The phase composition of the prepared samples was 
checked by x-ray diffraction (XRD) (Huber G670f cam- 
era, CuK Ql radiation, ImagePlate detector). The sam- 
ples contained BaCd(VO)(P04)2 and minor amount 
(~ 3%) of unreacted BaCdP20y. One then expects an 
equal amount of unreacted VO2 to be present in the 
samples. However, a minor impurity of VO2 can not 
be resolved by XRD due to the overlap of the strongest 
reflection of V0 2 with that of BaCdVO(P0 4 ) 2 . Instead, 
a minor amount of the VO2 impurity is evidenced by a 
small kink in the magnetic susceptibility data at 340 K 
(see Sec. Hn}. 

We tried to improve the quality of the samples by vary- 
ing temperature and duration of the annealing. Unfortu- 
nately, BaCdVO(P04)2 is rather unstable, and the for- 
mation of an unknown impurity phase was observed af- 
ter long annealings at 800 °C or any annealings above 
800 °C. The annealing at 900 °C resulted in melting and 
complete decomposition (as seen by powder XRD) of the 
compound towards unknown phases^ Regarding these 
difficulties, we did not attempt to grow single crystals 
of BaCdVO(P04)2, because the data measured on poly- 
crystalline samples are sufficient to determine the param- 
eters of the FSL model and the location of the system 
within the phase diagram. 

Magnetization (M) data were measured as a function 
of temperature using a SQUID magnetometer (Quantum 
Design MPMS). Specific heat C P (T) was measured on a 
pressed pellet with a standard relaxation technique us- 
ing a Quantum Design PPMS. All the measurements 
were carried out over a wide temperature range (0.4 
K < T < 400 K) and a field up to 7 T. The low- 



temperature measurements were done partly using an 
additional 3 He setup. 



III. RESULTS 

Figure [3] shows the magnetic susceptibility \ — M/H 
of BaCdVO(P04)2 as a function of temperature. At 
high temperature (above 20 K), the data behave in 
a Curie-Weiss (CW) manner. At lower temperature 
and low fields, the curve exhibits a broad maximum 
at T* ax ~ 2.7 K. The maximum is characteristic for 
low-dimensional spin systems and indicates a crossover 
to a state with antiferromagnetic correlations. A slight 
change of slope observed at ^ 1 K is associated with mag- 
netic ordering, while the Curie-like upturn in low fields 
below 1 K is likely caused by extrinsic paramagnetic im- 
purities. An increase of the magnetic field leads to a 
shift of the maximum to lower temperatures, an increase 
of the x value at the maximum (xmax), and an increase 
of x at the lowest investigated temperatures. This be- 
havior is quite similar to that observed in the other FSL 
systems 2 - and rather typical for low-dimensional and/or 
frustrated antiferromagnetic systems. 

Careful examination reveals a small kink in 1/x vs - T 
curve at 340 K. This kink is indicative of the VO2 im- 
purity, since VO2 undergoes metal- insulator transition 
at 340 KM- We emphasize that this impurity does not 
affect any of the results reported below, as the paramag- 
netic contribution of VO2 below 300 K is temperature- 
independent, while the energy scale of the exchange in- 
teractions in BaCdVO(P04)2 is lower by two orders of 
magnitude. However, the impurity prevents us from us- 
ing the data above 300 K in further analysis. 

High-temperature (20 — 300 K) susceptibility data can 
be fitted with a CW law corrected for the temperature- 
independent contribution xo that accounts for diamag- 
netism of core shells and Van Vleck paramagnetism: 



X(T) = Xo 



C 



T + 6>cw 



(1) 



The fitting resulted in xo = — 4-4(1) x 10 4 emu/mol, 
C = 0.375(2) emu K/mol, 6 CW = 0.66(4) K. This C 
value corresponds to an effective magnetic moment /x e ff = 
1.73(1) hb in perfect agreement with the expected spin- 
only value for V +4 . Note that #cw is slightly dependent 
on the lower limit of the data used due to the curvature 
in 1/x related to the maximum in xO^)- However, values 
of #cw below 1 K are obtained for any fit with a lower 
limit above 20 K. 

Thus, we find a very low 6qw value, despite the sus- 
ceptibility maximum at T^ ax ~ 2.7 K indicates some- 
what stronger exchange interactions. This result implies 
that several interactions with different signs are present 
in the system under investigation. To get a quantitative 
estimate of the exchange interactions, we fit the suscep- 
tibility data above 5 K with the high-temperature series 



4 



o 

! 

a 



25 



solution (a) 

solution (b) 

H fl = 0.5T 





FIG. 3: (Color online) Temperature dependence of the sus- 
ceptibility x(T) of BaCd(VO)(P0 4 )2 measured at different 
applied fields. Insets show inverse susceptibility (open cir- 
cles) along with the HTSE [Eq. (J2Jl, upper inset] and the 
Curie- Weiss [Eq. {l}, lower inset] fits. In the upper inset, 
dashed and solid lines denote solutions (a) and (b), respec- 
tively. In the lower inset, solid line is the Curie- Weiss fit, 
while the arrow indicates the kink at 340 K due to the VO2 
impurity. 



expansion (HTSE) for the FSL modeljH 



X(T) = Xo + 



N A g 2 ^ B 
knT 



E 



Ji 

knT 



(2) 

where Xo is temperature-independent contribution, g is 
Lande g-factor, and c m ^ n are the coefficients listed in Ta- 
ble I of Ref. Gj. We find xo = -3.8(2) x 1CT 4 emu/mol, 
g = 1.968(3), Ji = -3.62(5) K, and J 2 = 3.18(2) K [solu- 
tion (a)], or, alternatively, xo — —3.9(2) x 10 -4 emu/mol, 
g = 1.981(3), Ji = 2.16(1) K, and J 2 = -2.05(4) K [solu- 
tion (b)]j2i The fit with the solution (b) looks somewhat 
better and extends to lower temperatures as compared 
to the fit with the solution (a) (see the upper inset of 
Fig. [3]). Therefore, one may speculate that the solution 
(b) is the correct one. However, the difference between 
the two fits is caused by the lower exchange couplings 
in the solution (b), since the HTSE is valid at T > Ji. 
In the region used for the fitting (above 5 K), both the 
solutions produce the fits of similar quality, hence they 
are indistinguishable. The solution (a) locates the sys- 
tem in the CAF region close to the FM critical region 
(i.e., the region on the CAF-FM boundary), while the 
solution (b) corresponds to a position deep into the sta- 
ble NAF region. The presence of two solutions in fitting 
X(T) with the HTSE for the FSL is a well-known prob- 
lem. The ambiguity has to be resolved by the use of other 
experimental data. 

To discriminate the valid set of J's, we turn to mag- 
netization data and analyze the value of the saturation 
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FIG. 4: Magnetization (M) as a function of the applied field 
(H) measured at 0.5 K - well below Tjv. The arrow marks 
the saturation field fioH 3 . 



field H s . According to theoretical results by Schmidt et 
air the saturation field can be calculated as 

fi H s = BZ [ (l - I (cos Q x + cos Q y )) cos <p 



9VB 

-(1 — cos Q x cos Qy) sin tp\ 



(3) 



where z = 4 (magnetic coordination number), 5=1/2 
is the spin value, tp and J c are defined in Section [H and 
(Qx,Qy) is the wave vector of the ordered state. Using 
the appropriate wave vectors for the CAF and NAF re- 
gions, one finds (J,oH s = 2(Ji + 2J 2 )k B / (gfi B ) for CAF 
and AJiks / (gi^s) for NAF phases. The first set of J's 
[solution (a)] results in /j,qH s = 4.2 T, while the second 
set [solution (b)] gives rise to a much higher saturation 
field fi H s = 6.55 T. 

Figure 3] presents the magnetization curve measured 
at 0.5 K, i.e., well below the ordering temperature Tjv- 
At low fields, the M(H) dependence is linear, while a 
positive curvature is observed above 2 T leading to a 
marked kink at the saturation field ixqH s ~ 4.2 T. 35 The 
saturation magnetization M s ~ 0.9 /x_b/V +4 is slightly 
less than the expected value of 1 [Ib- The reduction of 
M s is likely caused by the weight error due to the pres- 
ence of non-magnetic impurity (BaCdP207) in the sam- 
ples under investigation (see Sec. [TT|'^ The curvature 
above 2 T is also somewhat puzzling. At first glance, one 
may interpret it as a spin-flop transition. However, the 
magnetic anisotropy of V +4 is weak. Thus, a spin-flop 
transition in Pb2VO(P04)2 is observed at 1 T^ and 
a strong increase of the anisotropy in the isostructural 
BaCdVO(P04)2 compound is unlikely. In the next sec- 
tion, we will suggest an alternative and plausible expla- 
nation for the curvature of the M(H) dependence above 
2 T. 

The experimentally determined saturation field 
matches exactly the value calculated for the solution 
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T(K) 

FIG. 5: (Color online) Temperature dependence of the specific 
heat C P (T) of BaCd(VO)(PC>4)2 measured at zero applied 
field. Open circles are the raw data, the solid line shows 
the phonon contribution [according to the fit with Eq. Q], 
and the dashed line indicates the magnetic contribution C mag . 
The anomaly at Tjv — 1 K is shown in detail in Fig. [6] In the 
inset, the magnetic entropy (S'mag) is plotted as a function of 
temperature, dashed lines show the entropy of 0.5/Zln2 and 
the respective temperature. 

(a) of the susceptibility fit but is far below the value 
expected for the solution (b). This clearly demonstrates 
that the solution (a), with a ferromagnetic J\ and 
an antiferromagnetic J 2 is the correct one. Thus, 
BaCdVO(P04)2 belongs to the frustrated ferromagnetic 
square lattice system as all the other AA VO(PC>4)2 
structural homologs. 

Specific heat (C p ) measurement at zero field is shown 
in Fig. [5] At high temperatures, C p is completely dom- 
inated by phonon excitations. Below 5 K, a decrease 
of temperature is accompanied by an increase of C p indi- 
cating that the magnetic contribution to the specific heat 
becomes prominent. At low temperatures, C' P (T) shows 
a broad maximum at T^ ax ~ 1.5 K due to the spin cor- 
relations and a small but well resolved peak at Tjv ~ 1 K 
associated with the magnetic ordering (see also Fig. [5]). 
Below Tjv, C p (T) drops rapidly. 

To extract the magnetic contribution to the specific 
heat C mag , we subtract an estimated phonon contribution 
from the total measured specific heat C p . For this pur- 
pose, the experimental data at high temperatures (15 K 
< T < 200 K) were fitted with 

(4) 

where A/T 2 accounts for the magnetic contribution at 
high temperatures, R — 8.314 J/mol K is the gas con- 



stant, and the sum of Debye functions accounts for the 
phonon contribution. The use of several Debye functions 

with distinct characteristic temperatures D is necessary 
due to the large difference in atomic masses of the ele- 
ments forming the BaCdVO(P04)2 compound. In gen- 
eral, the procedure is similar to that reported in Refs.l26l. 
[37l and [3a for the related vanadium phosphates. The re- 
sulting C mag (T) curve is shown in Fig. [5] In the inset, 
we plot the temperature dependence of the magnetic en- 
tropy S mag (T) as obtained by integrating C mag (T)/T. At 
high temperatures, S mag (T) converges towards the value 
Soo ~ 5.3 J/mol K, in reasonable agreement with the ex- 
pected value i?ln2 = 5.76 J/mol K taking into account 
the uncertainty in the estimated phonon contribution. 
The result indicates that C mag reflects the intrinsic con- 
tribution of BaCdVO(P0 4 ) 2 - 

At the lowest limit (T = 15 K) of the specific heat fit 
to Eq. (UJ), the magnetic contribution A/T 2 amounts to 
a small part of the total specific heat only. Therefore, 
the parameter A is obtained with a large uncertainty, 
A = 100 ± 30 J K/mol. The A/T 2 term corresponds to 
the lowest order in the HTSE for the specific healjii and 
may be expressed as A = 0.375J C = 0.375i?(J 1 2 + Jf ). 
We find A = 72 J K/mol and A = 29 J K/mol using the 
Ji values of the solutions (a) and (b) of the susceptibility 
fit, respectively. Both theoretical values are below the 
experimental one, but the result for the solution (a) is 
still within the error bar, while that for the solution (b) 
is far below. Thus, this analysis also favors the solution 
(a) and the CAF scenario for BaCdVO(P0 4 ) 2 . 

A further rough estimate of the thermodynamic en- 
ergy scale J c can be obtained from the T-dependence 
of the magnetic entropy. In a general approximation, 
J c is approximately twice the temperature at which the 
entropy reaches half of its high-temperature limit. In 
BaCdVO(P0 4 ) 2 , S(T) = 0.5i?ln2 at T = 2.5 K. This 
temperature corresponds to J c = 5 K, in perfect agree- 
ment with the solution (a) (J c = 4.8 K), but in poor 
compliance with the solution (b) ( J c = 3.0 K). 

Field dependent specific heat measurements are pre- 
sented in Fig. [5] The increase of the field results in the 
suppression of the broad maximum at 1.5 K, while the 
absolute value of C p at Xjy is enhanced. The maximum 
enhancement of the peak value is observed at 2 T where 
the transition anomaly is most pronounced. With fur- 
ther increase of the field, the peak value decreases again 
and the anomaly broadens slightly. However, between 3.5 
and 4 T the anomaly gets completely suppressed. The 
variation of the transition temperature is as follows: Jjv 
increases very slightly from to 2 T up to a maximum 
Tn(2 T) = 1.1 K and decreases more clearly with fur- 
ther increase of the field. This result enables to draw the 
H — T phase diagram of the system (inset of Fig. . The 
field dependence of the specific heat is similar to that 
of Pb2VO(P04)2r^ but the smaller energy scale of the 
exchange couplings facilitates experimental access to the 
larger part of the H vs. T phase diagram. 

The field dependence may be understood as follows. 
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FIG. 6: Specific heat of BaCdVO(P04)2 measured in different 
magnetic fields. The inset shows the H — T phase diagram of 
the system. 

The broad maximum at 1.5 K is caused by short-range 
antiferromagnetic correlations. These correlations are 
suppressed by the applied magnetic field, hence the max- 
imum is strongly reduced already at 1 T. The magnetic 
entropy is transferred to the transition anomaly resulting 
in a much larger and sharper peak at Tjv- As we have 
mentioned in the introduction, long-range magnetic or- 
dering in low-dimensional and/or frustrated spin systems 
is suppressed by quantum fluctuations. Magnetic field 
suppresses these fluctuations, therefore Tjv is slightly en- 
hanced at low fields. However, above 2 T the field is 
strong enough to overcome the antiferromagnetic order- 
ing, hence Tn is reduced and finally suppressed below 
0.4 K at 4 T. 



IV. DISCUSSION 

Our experimental results and our analysis demonstrate 
that BaCdVO(P04)2 is a frustrated square lattice with 
a ferromagnetic exchange J\ along the side of the square 
and an antiferromagmetic exchange J2 along the diago- 
nal. The first evidence for the presence of antiferromag- 
netic and ferromagnetic exchange of similar size is given 
by the very small value of the Weiss constant 6*cw < 1 K 
contrasting a much higher value of the temperature 
T^ ax ~ 2.7 K of the maximum in x(T'). The latter one 
reflects the onset of antiferromagnetic correlations. In 
a non- frustrated S = 1/2 antiferromagnetic square lat- 
tice, this maximum is expected at T* ax = 0.95J (Refs. 
l39l and , suggesting some antiferromagnetic exchange 
with J > 3 K in BaCdVO(P0 4 ) 2 . Then a much smaller 
$cWj which is the sum of all the exchange interactions 
in the system under investigation, implies an additional 
exchange of the opposite sign and thus ferromagnetic. 

A fit of the susceptibility in the T range 5 K < 



T < 300 K with the HTSE for the FSL resulted in 
two sets of exchange coupling parameters: solution (a) 
with Ji ~ —3.6 K, J2 ~ 3.2 K and solution (b) with 
J 1 ~ 2.2 K, J 2 — —2.1 K. The two solutions correspond 
to different regions of the FSL phase diagram and realize 
the frustrated CAF and non-frustrated NAF regimes, re- 
spectively. Further experimental results enable to select 
(a) as the correct solution. 

The frustrated regime of the solution (a) causes the 
low saturation field of hqH s = 4.2 T. The non-frustrated 
ground state for the solution (b) corresponds to the 
higher saturation field of /j,qH s = 6.55 T despite both 
J\ and Ji are weaker than in the solution (a). Magne- 
tization measurements at low temperatures show a well 
defined saturation at the saturation field \lqH s = 4.2 T. 
This value matches perfectly the noH s expected for the 
solution (a) but is well below the ^qH s expected for the 
solution (b). Further on, an analysis of the magnetic spe- 
cific heat at high temperatures (> 15 K) as well as of the 
T dependence of the magnetic entropy at low tempera- 
tures (< 3 K) leads to a thermodynamic energy scale of 
the exchange couplings J c in reasonable agreement with 
the solution (a) but in clear disagreement with the solu- 
tion (b). Therefore, the solution (a) is the appropriate 
one to describe the spin system of BaCdVO(P04)2- 

The solution (a) with J\ ~ —3.6 K and J2 — 3.2 K 
places the system to the CAF region of the FSL phase di- 
agram. With J2/J1 = —0.9 corresponding to ip/n = 0.76, 
BaCdVO(P0 4 ) 2 is quite close to the border of the FM 
critical regime ( J2/J1 = —0.7), closer than any of the pre- 
viously reported FSL compounds such as SrZnPO(P04)2 
(J2/J1 = -1.1) and Pb 2 VO(P0 4 ) 2 or BaZnVO(P0 4 ) 2 
{J2/J1 = —1-8). A direct evidence for the strong frus- 
tration in BaCdVO(P04)2 is given by the position of the 
specific heat maximum (Tj^ ax ) and the maximum value 
of the magnetic specific heat (C max ). The specific heat 
of the FSL at low temperatures (T < J c ) is not known 
precisely, since there is presently no established way to 
compute it. However, for an unfrustrated square lat- 
tice, the value at the maximum seems nowadays to be 
established within a few percent, C max ~ 0.46i? = 3.82 
J/mol K at T^JJx = 0.604M1 It is further established 
that tuning the FSL towards a critical region leads to 
a broadening of this maximum, a reduction of the C mag 
value at the maximum, and a shift of the maximum to- 
wards lower temperatures. Such a trend was already 
observed in the previously studied FSL systems, where 
C m a X /i? and T^ ax /J c are decreased from 0.46 and 0.6 in 
Li2VOSi04 (far away from any critical region) towards 
0.44 and 0.40 in Pb 2 VO(P0 4 ) 2 and finally to 0.40 and 
0.29 in SrZnVO(P0 4 ) 2 (closer to the critical region)^ 

In BaCdVO(P0 4 ) 2 , C max /i? reaches only 0.33, well 
below the value found in the other FSL systems, while 
T^/Jc = 0.31. The very small value of C max cannot be 
attributed to a scaling problem due to a large amount of a 
foreign phase, since the magnetic entropy at high temper- 
atures is close to R In 2 (the phonon contribution at T^ ay . 
amounts to less than 1% of the total specific heat and can 
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therefore safely be neglected). Further on, a comparison 
of C* ma g(T) near its maximum for different FSL systems 
using a reduced T/J c temperature scale indicates that 
the maximum in C mag (T) of BaCdVO(P04) 2 is broader, 
and the decrease of C mag (T) above T^ ax is much weaker 
than in the other known FSL systems. Thus, this broad 
maximum with a small C mag value at the maximum is 
also a direct experimental evidence for the strong frus- 
tration. 

A further (and new) manifestation of the frustration 
is likely observed in our magnetization curve (Fig. [?]). 
Above 2 T, the slope of M(H) increases with magnetic 
field and steeps up just before saturation is reached. Such 
a behavior has been recently predicted by theoretical cal- 
culations for the FSLfii This steeping up just below H s 
is suggested to be more pronounced for systems that are 
close to the critical regions. However, further magnetiza- 
tion measurements on other FSL systems are needed to 
confirm this trend. 

Using exact diagonalization of finite-size clusters, 
Shannon et calculated the dependence of a number 
of characteristic properties of the FSL (such as Tj* ax , 
7^, Xmax, and C max ) as a function of ip. Compar- 
ison of these predictions with the experimental obser- 
vation for BaCdVO(PC"4)2 and the related compounds 
indicates that while the overall if dependence seems to 
be reproduced, the calculated absolute values are slightly 
too large. This difference is likely a consequence of the fi- 
nite cluster size that leads to a gapped energy excitation 
spectra and thus an exponential vanishing of C mag (T) 
and x(T) at very low T in contrast to the real behavior. 
The missing entropy has then to be recovered at higher 
temperatures, thus enhancing C mag (T) at its maximum 
and shifting the maximum to higher T. 

A last comment on the antiferromagnetic order ob- 
served at TV — 1.0 K. The anomalies in x{T) and C P (T) 



are much smaller than those expected for a classical 
three-dimensional magnetic system but very similar to 
those observed in the other FSL systems, e.g.JLi^VOSiCU 
(Refs. and or Pb 2 VO(P0 4 ) 2 (Refe. M and [H). 
Since in some of these systems long-range magnetic or- 
der was directly confirmed by NMR or neutron scattering 
experiments ; 18 ' 2 ^ 28 the similarity in the behavior allows 
us to safely claim that these small anomalies in x(T) 
and C P {T) correspond to the onset of (columnar) anti- 
ferromagnetic order in BaCdVO(P04)2 too. The ratio 
R = Tn/T*^ is smaller in BaCdVO(P0 4 ) 2 (R = 0.39) 
than in the other FSL systems (0.4 < R < 0.6), reflect- 
ing either the weakness of intcrlayer exchange and/or the 
suppression of the AFM order by quantum fluctuations. 

In conclusion, our study of BaCdVO(P04) 2 shows that 
the magnetic properties of this compound are well un- 
derstood within the FSL model. Magnetic susceptibility, 
magnetization, and specific heat data consistently sug- 
gest J\ ~ —3.6 K, J 2 ~ 3.2 K, hence a ~ —0.9, lo- 
cating the compound into the CAF region of the FSL 
phase diagram. BaCdVO(P04) 2 lies closer to the criti- 
cal region of the FSL than any of the previously reported 
compounds. This conclusion is supported by a strongly 
reduced maximum of the magnetic specific heat and a 
positive curvature of the magnetization curve consistent 
with the recent theoretical predictions. 
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